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Linear perturbation theory is appropriate to describe small oscillations of stars, while a mild 
non-linearity is still tractable perturbatively but requires to consider mode coupling, i.e. to take 
into account second order effects. It is natural to start to look at this problem by considering the 
coupling between linear radial and non-radial modes. A radial pulsation may be thought of as an 
important component of an overall mildly non-linear oscillation, e.g. of a proto-neutron star. Radial 
pulsations of a spherical compact objects do not per se emit gravitational waves but, if the coupling 
between the existing first order radial and non-radial modes is efficient in driving and possibly 
amplifying the non-radial oscillations, one may expect the appearance of non-linear harmonics, and 
gravitational radiation could then be produced to a significant level. More in general, mode coupling 
typically leads to an interesting phenomenology, thus it is worth investigating it in the context of 
star perturbations. 

In this paper we develop the relativistic formalism to study the coupling of radial and non-radial 
first order perturbations of a compact spherical star. From a mathematical point of view, it is 
convenient to treat the two sets of perturbations as separately parametrized, using a 2-parameter 
perturbative expansion of the metric, the energy-momentum tensor and Einstein equations in which 
A is associated with the radial modes, e with the non-radial perturbations, and the Ae terms describe 
the coupling. This approach provides a well-defined framework to consider the gauge dependence 
of perturbations, allowing us to use e order gauge-invariant non-radial variables on the static back- 
ground and to define new second order Ae gauge-invariant variables representing the result of the 
non-linear coupling. We present the evolution and constraint equations for our variables outlining 
the setup for numerical computations, and briefly discuss the surface boundary conditions in terms 
of the second order Ae Lagrangian pressure perturbation. 

PACS numbers: 04.30.Db, 04.40.Dg, 95.30.Sf, 97.10.Sj 



I. INTRODUCTION 

Since the discovery of pulsars in the sixties neutron stars have acquired a special status in physics: as supernovae 
remnants they are fundamental to our understanding of the final stages of evolution and fate of upper main sequence 
stars; as the most compact observed objects they are a test-bed for strong-field gravity, i.e. general relativity and 
its generalizations; they are a unique laboratory for fundamental physics such as nuclear interactions, superfluidity, 
superconductivity. Either as isolated objects or in binary systems neutron stars are important gravitational wave 
sources, and in the near future the analysis of gravitational radiation from neutron stars will open up a direct window 
on their interior, possibly revealing details on the equation of state of nuclear matter, the dynamics of the crust-mantle 
interaction and the inner superfluid/superconducting core. Early stages of neutron star formation during iron core 
collapse in supernovae and other possible transient dynamical phases are particularly interesting from the point of 
view of gravitational wave physics. In this case the gravitational radiation emitted should carry the signature of the 
yet not understood mechanism driving the explosion and that of the unknown phenomena responsible for the observed 
kick velocity of many pulsars, providing complementary information to that carried by neutrinos. 

Physics of compact objects like supernovae core collapse or neutron stars may be studied with various approaches 
and approximations, from purely Newtonian, or Newtonian with relativistic corrections, to the full general relativistic 
treatment of numerical relativity. For many astrophysical problems even a purely Newtonian approach may be 
adequate. However, in order to study neutron stars as sources of gravitational waves, a full general relativistic 
treatment is more satisfactory because in this case the gravitational radiation is built-in in the calculations rather 
than calculated a-posteriori with the quadrupole formula, thus one can obtain more accurate results. As for other 
gravitational wave sources, in the long term the goal is a full numerical relativity treatment coupled with a detailed 
description of the matter physics. Meanwhile, relativistic perturbation methods remain a valid alternative: they are 
computationally far less intensive because, even in dealing with non-linear effects as we shall be doing here, one has to 
solve linear partial differential equations instead of the complete non-linear Einstein equations. Hence, accurate results 
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can be obtained at a relatively low cost. Moreover, even when numerical solutions of the complete field equations 
are available, analytic perturbative methods often help to shed light on the physical processes at work, clarifying the 
interpretation of the numerical results. 

Linear perturbations and instabilities of neutron stars have been studied for long time 0, 0] but relatively little 
is known of non- linear dynamical effects (see e.g. 0, 0, 0j S 0) and therefore second order studies may help to 
understand known problems and even reveal a new phenomenology. Non-linear effects are the rule rather than the 
exception in phenomena in all branches of physics, and have to be taken into account for accurate modeling. This 
should even be more natural when the mathematical modeling is not phenomenological, but rooted in a fundamental 
theory which is per se non-linear, as general relativity is. When we use perturbation theory, we generally expect that 
if the perturbations are very small, then second order effects should be negligible. On the other hand, while in order 
to treat strong non-linear dynamical effects a fully non-perturbative approach is required, one may expect that much 
of the interesting physics only involves a mild non-linearity for which a second order treatment should be perfectly 
suited. Furthermore, a proper parametrization of a problem may even lead to the unexpected result that a second 
order treatment appears to provide a good description of the physics even in a mildly non-linear regime where a priori 
one would expect the perturbative approach to fail. An example is given by the study of black hole collision in the 
close limit approximation 8]. With this in mind, it is reasonable to expect that a second order treatment of non- linear 
oscillations of neutron stars should be adequate in many circumstances of astrophysical and/or gravitational interest. 

In this paper we develop the formalism to study a specific non-linear effect, focusing our investigation on the 
coupling of radial and non-radial first order perturbations of a relativistic star. Work is currently under way to 
apply this formalism and carry out the actual study, and will be presented elsewhere 0| . The main goal is that of 
understanding whether this coupling can lead to new effects such as fluid instabilities and/or persistence of oscillations 
and amplifications that can produce a significant amount of gravitational radiation. Furthermore, this work may be 
a first step toward a more comprehensive study of second order perturbations of compact stars and mode coupling. 

The physical picture we have in mind here is that of a star undergoing a phase of overall oscillation or wobbling 
that we want to describe going beyond the linear approximation. Purely radial modes are going to be a natural 
component of such dynamical phase; on the other hand, even leaving rotational effects aside (as we do here) it is hard 
to think of oscillations that don't have at least a tiny non-radial component. While purely radial and non-radial second 
order effects may also eventually become important in the non-linear phase, it is thus natural to first investigate the 
coupling of their linear components. A specific motivation for this is that the purely radial oscillations of a spherical 
star don't emit per se any gravitational waves but, if the coupling we aim to study is efficient in driving and possibly 
amplifying the non-spherical oscillations, one may expect the appearance of non-linear harmonics, and gravitational 
radiation could then be produced to a significant level, provided that at first order both the radial and the non-radial 
modes are non- vanishing. More in general, mode coupling typically leads to an interesting phenomenology, thus 
worth investigating. For example, a specific effect we can easily anticipate is that axial modes, decoupled from fluid 
perturbations at first order, can be sourced by radial oscillations. Finally, a further motivation is that there are a 
number of studies aiming at investigating if non-radial oscillations of stars can be excited by external sources (see 
6-g- Instead, our general idea is to see if the non-radial oscillations can be driven or even amplified through 

coupling by an internal radial oscillation, regardless of the presence of an external source. These sort of non-linear 
processes could occur, for instance, in a proto-neutron star that is still pulsating. A mainly radial pulsation could for 
example drive the non-radial oscillations, either naturally present, or excited through fall-back accretion. 

Our study can also be seen from a slightly different perspective: we are investigating linear non-radial perturbations 
of a radially oscillating star, where the radial oscillation is also treated perturbatively. From this point of view our 
problem is mathematically similar to that of st udy ing linear perturbations of a slowly rotating star (see e.g. 13j) 
described by the Hartle-Thorne metric [13 . ITsl HQj : in the latter case the spacetime we perturb is itself a special 
perturbation - stationary and axisymmetric (typically up to second order) - of a spherical static star. Similarly, from 
this point of view we consider here perturbations of a spacetime which is itself a spherically symmetric time-dependent 
perturbation of a static spherical star. As we shall see, this point of view turns out to be of practical value, although 
it must be used with care because of the typical gauge issues of relativistic perturbation theory. 

Mathematically is however more satisfactory to consider radial and non-radial perturbations as different first order 
perturbations of a static spherically symmetric background that couple at second order. This point of view renders 
transparent two crucial aspects of our problem: i) the perturbations are defined as fields on a static spherical back- 
ground; ii) the two sets of perturbations are separately parametrized. Thus a well-defined framework for our study is 
provided by a multi-parameter relativistic perturbation formalism previously developed 0, |0| . This allows us to set 
up the formalism in a hopefully transparent way, properly considering the gauge dependence of perturbations. Fixing 
the gauge for radial perturbations , we borrow the formalism developed by Gundlach and Martin Garcia |l9ll2(i [ (based 
on that of Gerlach and Sengupta |22; we shall refer to the GSGM formalism in the following), which gives equations 
for gauge-invariant perturbations on a general time-dependent spherical background. This then allows us to: i) have 
gauge-invariant non-radial first order variables on our static background and ii) to define new second order variables, 
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describing the non-linear coupling of the the radial and non-radial linear perturbations, that are also gauge-invariant 
at second order. This higher order gauge invariance, attained by partially fixing the gauge at first order, is similar to 
that considered for example in and ^ , although in our case we deal with a 2-parameter expansion 0, and 
we only need to fix the gauge for radial perturbations. 

At first order most of the fluid perturbations appear in the polar part of the spectrum. Hence we expect that the 
effects of the radial non-radial coupling will predominantly manifest themselves through polar modes. Therefore in 
this paper we focus on polar perturbations of a perfect fluid star. 

The plan of the paper is the following. In Section II we describe the general 2-parameter perturbative framework 
we are going to use. In Section III we briefly recall the GSGM formalism. Using this later, in Section IV we introduce 
the radial and non-radial first order perturbations of a static, spherically symmetric star. Section V introduces the 
second order perturbations that account for the coupling between the radial and non-radial first order ones; we prove 
the gauge invariance of such perturbations, and give the equations they fuUfiU. In Section VI we briefly discuss the 
problem of the boundary conditions, defining the second order Lagrangian pressure perturbations needed to fix the 
problem at the surface of the star. Finally in Section VII we draw our conclusions. 

The conventions that we follow throughout this work are: Greek letters are used to denote spacetime indices; 
capital Latin letters are used for indices in the time-radial part of the metric; lower-case Latin indices are used for 
the spherical sector of the metric. We use physical units in which SttG = c ~ I. 



II. PERTURBATIVE FRAMEWORK 



In order to study the coupling of the first order radial and non-radial perturbations of a static spherically symmetric 
star it is convenient to use a 2-parameter perturbative approach, where the coupling will then appear at higher order. 
A natural framework for the study of this problem is therefore provided by the multi-parameter non- linear perturbative 
formalism introduced in |l7l llSll. Generalising well known mathematical ideas at the basis of standard 1-parameter 
linear (23. 12^ and non-linear |26l |2^ perturbation theory, the basic underlying assumption in the construction of a 
multi-parameter relativistic non-linear formalism is the existence of a multi-parameter family of spacetime models 
that can be Taylor expanded around a background spacetime, representing an idealized situation. These spacetime 
models are labeled by a set of parameters that formally control the strength of the perturbations with respect to the 
background, and serve as book keeping. The crucial point to obtain a manageable theory is to choose a convenient 
background, in our case the static spherically symmetric star. 

Let us denote the metric of this background with g(*''0), i.e. a Tolman-Oppenheimer-Volkov (TOV) solution of the 
field equations. Denoting with g the physical metric, we shall expand it in the two parameters A and e, and we will 
use superscript indices j) to denote perturbations of order i in A and j in e. 

We then have 
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where the terms g(i'°) and g(°'i) respectively represent first order radial and non-radial perturbations, and g(i'^) is the 
non-linear contribution due to the coupling, which is the new quantity that we want to compute. The other second 
order perturbations, i.e. the self-coupling terms of order A^ and e^, will not be considered in this work. Any field can 
be expanded as the metric in Eq. JQ). In particular, we can expand in this way fluid variables like the energy density 
and the 4-velocity, and for the energy-momentum tensor T we can formally write 



(2) 



where each term T^^'-'^ collects terms of the metric and fluid variables of the appropriate order. Let us now consider 
the structure of the perturbed field equations, following a standard procedure ,2j|. We start from thcfuU Einstein 
equations: 



£;[g,^^] = G[g]-T[g,t/;^] = 0, 



(3) 



where G denotes the Einstein tensor, and i/?^ (A=l, . . .) the various fluid variables. If we introduce the perturbative 
expressions and (0) into Eq. we can expand the latter up to (1, 1) order, obtaining 



.(0,0)^^(0,0) 



;^^(l,0) 
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(1,0)^^(1,0) 



r(0,l) 



(0,1) 
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The previous equation is satisfied for arbitrary values of the two parameters if, and only if, each coefficient of the 
expansion vanishes. Therefore, setting each of these terms to zero, E^'^'^'' = represents the TOV equations (see, 
e.g., 28J), while each of the other = terms represent the perturbative equations of order (z, j). As differential 

operators, the act linearly on each of the terms in square brackets, while they are non-linear functions of the 

background quantities g^"'") and ■0^'*''' . 

At first order in A, we obtain the equations describing the radial perturbations on the TOV background, 



.(1,0) 



.(1,0)^^(1.0) 



= 0. 



(5) 



Linear radial perturbations have been extensively analyzed in the literature (see |28l | and references therein, and |2£ 
for more recent results). 

The linearized equations for the non-radial perturbations come from the first order terms in e, 



^(0,1) g(o-i), ^(^0.1)1 - 



= 0, 



(6) 



and were first studied by Thorne and Campolattaro (Sfl- Later they became the subject of many investigations, see 

e.g. pm HI a Hi III. 

Finally, the equations describing the radial non-radial coupling, the ones we shall focus on, have the following form 



£;(i,i) 



.(M),^(M) 



.(1,0) 



(0.1)^^(1,0) 



.(0,1) 



0. 
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It is an intrinsic feature of perturbation theory that the procedure to solve the above equations is iterative. Thus, when 
we arrive at the stage of solving Eq. Q, the terms g(i'°) and V'a ^i'^ assumed to be known from solving the radial 
linear equation while gC'^) and f/'^'^'' are solutions of the non-radial linear perturbation equations ©. Hence, 
because of the linearity of the operator E^^'^^ in acting on each of the terms in square brackets, these terms play the 
role of sources in Eq. Q . Taking into account the nature of the different sets of perturbations we are considering, it 
turns out that the operator E^^'^^ acts on the pair g'-^'^\ "^a'^^ ^'i- O same way that E^'^'^'' acts on g*-"'^-*, 

■j/)^'^-* in Eq. ©. The reason for this is that both operators E^'~''^^ and E^^'^'' come from the linearization, around 
the static background, of the Einstein tensor operator acting on non-radial perturbations. Therefore, using again the 
linearity of e'--^'^\ we can define 



as the non-radial perturbation operator. Hence, we can re-write equation |(7J) in the final form 



- r(o,i) 



g 



(1,1) 



.(1,0) 



.(0.1) 



.(0.1) 



(8) 



(9) 



The particular structure of these equations is very helpful in order to develop a numerical code for studying the 
coupling of radial and non-radial perturbations. We want to emphasize here that our goal is that of solving the 
perturbative equations in the time domain. To this end it is then very useful to rely on well known initial value 
formulations for the linear non-radial perturbations. Thus, assuming we have working numerical codes for the first 
order perturbations in A and e, in order to implement a code for the coupling Ae variables we only need to modify the 
code for the e variables by adding the source terms to the right-hand side of the evolution algorithm. Then, at every 
time step in the evolution, and having fixed the background, we have to: i) evolve the equations to obtain the value 
of the first order radial and non-radial perturbations; ii) to use their values to evolve the coupling variables. 



III. SUMMARY OF THE GSGM FORMALISM 



In this Section we briefiy recall the GSGM formalism, introduced by Gerlach and Sengupta 0, [s^l and further 
developed by Gundlach and Martfn-Garci'a 0, [2p, "SS] , to study first order gauge-invariant perturbations of a general 
time-dependent spherically symmetric stellar background. From the formal point of view of parameter expansion of 
Section m here we are dealing with standard 1-parameter linear perturbations of the form 



?./3 = gS3+egS, (10) 



where g^^^ is the time-dependent spherically symmetric background and the e perturbations are non-radial. 
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A. The time dependent perfect fluid background 

The background manifold is the warped product x S^, where S'^ denotes the 2-sphere and a two-dimensional 

Lorcntzian manifold. The metric can be written as the semidirect product of a general Lorentzian metric on M^, 
gABi and the unit curvature metric on 5*^, that we call ^ab'- 



^ 9ab 



(11) 



Hereafter, denotes the coordinates on and the coordinates on 5^; r = r(x^) is a function on that 
coincides with the invariantly defined radial (area) coordinate of spherically-symmetric spacetimes. A vertical bar 
is used to denote the covariant derivative on Af^ and a semicolon to denote the one on 5^, thus we have gAB\c = 
lab:c = 0. One can introduce the completely antisymmetric covariant unit tensors on and on 5*^, tAB and tab 
respectively, in such a way they satisfy: eAB\c = ^ab-.c — , ^ac^^^ = —9a ' ^^'^ ^ac^^'^ = ~7a • 

In this paper we consider a perfect-fluid description of the stellar matter, thus the energy-momentum tensor is 

= {P + P)UaUl3 + Pga0 ' (12) 

where p and p are the energy density and pressure, and Ua is the fluid velocity. 

In the spherically symmetry background taj3 has the same block diagonal structure than the metric, 

ta0 = diag ( tAB , r'^Q{x'^)lab) , (13) 

and the fluid velocity takes the form Ua = {ua, 0). An orthonormal frame on the submanifold can be formed from 
and the spacelike vector 

nA = -cabu^ => uau^ = . (14) 
The metric g^^ and cab can be written in terms of these frame vectors as follows 

Qab = -uaub + riAriB , sab = haub - uaub ■ (15) 

Then we have 

tAB = pUAUB + priAriB , Q=P- (16) 
In any given coordinate system for M^, {x^} , one can define the following quantity: 

VA = -r^A ■ (17) 
r 

Then, any covariant derivative on the spacetime can be written in terms of the covariant derivatives on and 5^, 
plus some terms due to the warp factor r^, which can be written in terms of va- Finally, the frame derivatives of a 
generic scalar function / are defined by 

f = u^f\A, f = n^f\A, (18) 
and we introduce the following background scalars: 

O = Inp, U = u-^VA, W = n^VA, M = u^\a^ v = n^^^ . (19) 

B. Perturbations 

Linear perturbations of a spherically-symmetric background can be decomposed in scalar, vector and tensor spherical 
harmonics. The scalar spherical harmonics F'"* are eigenfunctions of the covariant Laplacian on the sphere: 

^abyhn ^ ^ ^)yim _ (20) 

A basis of vector spherical harmonics (defined for 1) is 

-w-lm — \rlm aim — ^bxrlva /Ol ^ 
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where the F^^'s have polar parity (they transform as (—1)', hke the scalar harmonics, under parity transforma- 
tions, and are also called even-parity type) and the S'^""s have axial parity (they transform as (—1)'"'""'^ under parity 
transformations, and are also called odd-parity type). A basis of tensor spherical harmonics (defined for ^ > 2) is 

\rlra V^^^, rylra -v^/m i ^(^ ^ \^lra„, qItu nlra i qItu /OOA 

-'ab = lab, ^ ab — ^ -.ab " ^ lab. Jab = ^a-.b + ^b-.a ^ l^^7 

where the yj™ , Zj^ have polar parity and the S*^™ have axial parity. In this paper we will only consider perturbations 
with polar parity. 

The perturbations of the covariant metric and energy-momentum tensors can be expanded in this basis as 

h,lrn \^lm l-,lm \^lrn 

''■AB ^ "A 

5go.fi =1 I , (23) 

6taf3 =1 I . (24) 

- :ab 



Let X be an arbitrary tensor field on the background spacetime and dX its linear perturbation. It is well-known that 
under a a first order gauge transformation, generated by a vector field ^ living on the background, the perturbation 
of X transforms as 

SX ^SX + C^X . (25) 

Then, the perturbation SX is gauge-invariant if and only if the Lie derivative of the corresponding background 
quantity X with respect to an arbitrary vector field f vanishes: C^X = [s^. Using this well-known result, it is 
possible to show that a complete set of gauge-invariant variables, combinations of the perturbations Hab, ^a, K, G, 
6tAB, St A, 5t^ , 5t^ , is given by the following quantities: 

kAB = hAB - {PA\B +Pb\a) , (26) 
k = K-iv'^pA, (27) 
Tab = StAB ~tAB\cP° -tAcpfB-^BcpfA^ (28) 

t3 = ,5^3 - p^{Q\c + 2Qvc) + ^^^-^ Q G , (29) 

Ta = StA-tAcp"^ ~'^QG\A, (30) 
= 5f -r^QG, (31) 
where Ta is defined for ^ > 1 , is defined for Z > 2, and 

PA = hA- \r^G\A ■ (32) 

Therefore, any linear perturbation of the spherically-symmetric background can be written as a linear combination 
of these gauge- invariant quantities. Is is important to stress the fact that although we are considering a perfect fluid 
energy-momentum tensor, the gauge invariant quantities introduced above are defined for any energy-momentum 
tensor. The equation of state for a perfect fluid has the formp = p{p, s), s being the specific entropy. The corresponding 
sound speed, Cs, can then be defined through the thermodynamical derivative 

Throughout this paper we will consider the particular case of a barotropic fluid (i.e. constant specific entropy), hence 
we will have p = p{p) . 

The polar fluid perturbations are described as follows. The perturbations of the fluid velocity components can be 
written as 



Sua — 



7 nA + ^hABU 



ylra ^ ^Iraylra j ^ (34) 
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where a is defined for Z > 1 . The energy density and pressure perturbations can be cast in the following form (using 
the barotropic equation of state) 

Sp^CjpY^"", Sp^clSp. (35) 
In terms of these quantities it is possible to define a gauge-invariant set of fluid perturbations: 

a — a — p^ub , (36) 



A 

7 = 7 — n 



P^ua\b + ^u^{Pb\a ~ Pa\b) 



(37) 



U! = UJ ■ 



P^i^\A ■ (38) 



The tensor kAs can be decomposed in the frame {u^, n^}: 

kAB = rj{-UAUB + nAriB) + 4'{uaub + uaub) + tpiuAriB + tiaub) , (39) 

where 77, (f) and ip are scalars. It is useful to consider the following new scalar variable 

X = - fc + 77 , (40) 

in the place of (j). Then, combining Einstein equations with the energy-momentum equations we can obtain the 
following set of equations: for Z > 2, 



for I > 1, 



77 = 0, (41) 



(42) 
(43) 
(44) 
(45) 
(46) 



+ 2(/i - U)ij' 




Ik" - 2clUyj' 


= Sk , 




= , 


167r(p -I- p)a 




—a 


— Set 1 








S-y . 



And finally, for I > 0, 



(47) 
(48) 



87r{p + p)j = [ky + C^, (49) 
Snpuj = -k" + 2C/Vy + Cu,, . (50) 

where the expressions of S^, S^, Ca, Sa, S^^, S^, C^, can be found in [T^ . 

IV. RADIAL AND NON-RADIAL PERTURBATIONS OF STATIC RELATIVISTIC STARS 

In what follows we summarize the first order perturbative analysis of the oscillations of a perfect-fluid static star. 
As anticipated in Section ^ we consider separately the radial pulsations, parametrized by A, and the non-radial 
oscillations, parametrized by e. Therefore, consistently with the notation of Section m we will now explicitly use the 
indices («,j); however, to simplify the notation, from now on we will use a bar to denote quantities associated with 

the static spacetime, the background of our 2-parameter perturbative formalism. Thus we have g^^ = g^"^'^'', and in 
the same way p, p, and Ua- 

The equilibrium configuration is described by the static spherically-symmetric metric: 



=Q/3 
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with a perfect-fluid energy-momentum tensor 

tap = {p+P)UaUl3+Pga0, (52) 

where = (-e*, 0,0,0). The mass function is introduced by means of the equahty e'^^^'') = 1 - 2M{r)/r . Then, 
the TOV equations are: 

_ M + 4TTpr^ _ p^r 
^■^ ^ r{r- 2M) " JTp ' ^ ^ 

iV/r 47rpr2 . (54) 

Specifying the equation of state of the stellar matter one obtains a 1-parameter family of solutions of these equations, 
depending on the central density. For the barotropic equation of state we use, p — p{p), the background sound speed 
is = dp/ dp . 



A. Non-radial perturbations in the GSGM formalism 

The equations for the first order non-radial perturbations have been known for a long time [sol Isil Is^ . In [s^ they 
have been presented in the framework of the GSGM formalism. We do not write their explicit expressions here, because 
they can be obtained as a particular case of the equations that we will write in Section FV GI for the coupling perturbative 
terms of order (1, 1), like g^"'"'"'"-', considering the homogeneous part of those equations, i.e. simply neglecting the source 
terms. One can also obtain the same equations directly from the general GSGM equations H41I50|I . by considering 
the special case of a static background, represented by the following quantities, 

S^ = (e-*,0), n^ = (0,e-^), (55) 

A2 = f/ = 0, i> = $', , (56) 

r 

and where the frame derivatives of a scalar function / take the special form /' = e~^f^r and / = e~'^f,t. Boundary 
conditions have to be imposed at infinity, at the stellar origin and on the stellar surface. They can be found in the 
references cited above. 



B. Radial perturbations in the GSGM formalism 



Radial perturbations can be seen as a particular subcase of non-radial ones, namely those corresponding to the 
harmonic I = 0. The fact that the GSGM quantities are not gauge-invariant in the I = case does not represent a 
problem for our study since, when we will consider the gauge invariance of our (1, 1) variables in Section IVB I we will 
assume that the gauge of the perturbative order (1,0) has been somehow fixed. In practice, we fix the radial gauge 
as in pLOf. 

^(1.0)^0, A;(i-°^=0, (57) 

which considerably simplifies the equations. In this gauge, the metric perturbations have a diagonal form as in the 
background: 

gir=diag(/.i^i,°\0), (58) 



h'AB = v'^''°hAB + ^'^'^°\uAUB + nAns) . (59) 



Then, the components of /i^^'' are given by 



h 



(1.0) _ 
AB — 



diag (e2* (x(i'°) - 27y(i'°)) , , (60) 



where x^^'°'' = Tyt^'O) + (/)(i'0). 
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For the fluid velocity, from Eq. (|34() . we have 



Su 



(1:0) 



The other fluid perturbations are given by: 



(61) 
(62) 



C. Equations for the radial perturbations 

The equations for the quantities t^;(1'0), 77(1^0) ^^(1,0) ,^(1,0) describing the radial perturbations can be found in [T9l |. 
For our purposes it is more convenient to use a different set of variables. First, instead of using oj^^''^) we use the 
enthalpy perturbation H'-^'^\ which significantly simplifies the equations. The second change consists in replacing 
the metric perturbation x''^'*'' with the quantity S^^'^\ in order to use a set of variables consistent with the one we 
will use for the non-radial perturbations in Section IV CI The definitions of ij(i'°) and S'^^'^^ are: 



Him 



p + p P + P ' r 



(63) 



Furthermore, using equation (34) in we end up with a set of evolution equations for S'^^'"-', iJ^^'^) and ^(i''^) which 
does not contain the quantity r]i'^'^\ 



Him 



.Him- 



in 
W 



(1:0) 



Sim = ^8n{p + p)^.(''°\ 
rW 



and with the Hamiltonian constraint 



1 

r I cj 

The quantity //(i^'^) can be found, in terms of (S'^^'"' , ij'^^'^), 7(1-°)), from the following equation: 



M/jy^i'"'' = 47r(p + p) 



■5(1^°)+ 1 + 



1 



/y(i,o) 



The system of equations ()64I67|) is equivalent to the one used by Ruoff ,4Q] , cf . also [28 



(64) 
(65) 
(66) 

(67) 
(68) 



D. Boundary Conditions for radial perturbations 

Boundary conditions must be fixed at the center and at the stellar surface r — R. The latter is given by the 
vanishing of the Lagrangian pressure perturbation, Ap = 0. Following we can express Ap*^i'°) in terms of the 
radial renormalized displacement function (: 



Using the relation between C and 7(1^°) 



r^Ap^'-'^=-{p + p)cle-^^ 



we arrive at the following boundary condition on the surface: 



= 0. 



(69) 
(70) 
(71) 



r=R 



The behaviour of Si^'°^ and Hi^'°^ on the surface can be found from the general evolution equations (|66|l and H64|l . 
The analysis of the regularity conditions at the origin (r — 0) leads to the following expressions: 

(^(1,0) ^ ^(1,0) ^ ^(1,0) ^ ^(1,0) ) _^ (^(1,0) ^ + 0(^3 ) ^ ^(1,0) + oir^), H^m (t) + 0{r^), 7^^'°^ {t)r + 0{r')). (72) 
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V. COUPLING OF RADIAL AND NON-RADIAL PERTURBATIONS OF RELATIVISTIC STARS 

At this point we have already specified the equations that determine the static background metric g(*'^'^), the radial 
perturbations of order A, and the non-radial perturbations of order e. The next crucial step is to find the equations for 
the perturbative coupling terms of order Ae, like the metric perturbations g(^'^). To this end we can further expand 
the first order 1-parameter GSGM formalism of Section IIIII as follows. The time dependent spherically symmetric 
variables can be split into a static and a radially oscillating parts, e.g. for the metric 

JO) _ (0,0) . (1,0) 

The non-radial first order perturbations on this time dependent spacetime can be split into a part which is the non- 
radial perturbation of the static background and a further term that describes the coupling. For the metric we then 
have 

Inserting (|73|l and H74|) into Eq. H10|l we re-obtain the 2-parameter expansion , and similarly for the fluid variables. 

This approach, based on splitting the GSGM variables and their equations, is very convenient because it saves 
us from the need of computing directly a second order expansion of Einstein's equations. In addition, following this 
approach we shall show in Section fV BI how to build gauge-invariant quantities associated with the perturbations g^^'-'^^ 
and the fluid variables of order (1, 1). 



A. GSGM formalism on a radially oscillating star 



In order to implement the GSGM formalism on a radially oscillating star, which is itself treated perturbatively, we 
split the quantities u^, n^, U, W, /i, in a static background part and a radial perturbation, as illustrated above. 
The frame vector fields are given by 



(1.0) 



.Ae-V'") 



(Ae-V-,(l-A4^)e-) 



The frame derivatives of a scalar perturbation, f^^'' = Z^"'^-* + A/*^^'^-*, on the radially oscillating star are: 

fW' ^ ^ g-A ^(0,1) ^ ^ |g-A ^ g-*^(l,0) ^(0,1) _ g-A ^(^0,1) 

The remaining quantities describing the spherical star are 



e ^ (1,0) 



U = u'^VA = A 7V 

r 



(1.0) 

W = n'^VA = ( - - A-" 

r 2r 



1 ^X 



(1.0) _ ^^(1,0) 



ll.r 



1 



$_^^(1,0) 



(75) 
(76) 

(77) 
(78) 

(79) 
(80) 
(81) 
(82) 



where in 1)81(1 we have used the ((66|l . 



B. Gauge invariance of the Ae coupling perturbations 



In this Section we are going to construct a set of gauge-invariant quantities at the (1, 1) perturbative order. As we 
shall see, we can find them with the help of the GSGM formalism. The main idea behind our construction is to build 
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the (1, 1) gauge-invariant variables starting from the gauge-invariant quantities of the GSGM formahsm, considering 
these latter as perturbations of a radially pulsating spacetime which is itself described as a perturbation of the static 
background, i.e. gf^ = g^p^'' + ^Sa^^^ ■ More precisely, let C/^^^ be any of the gauge-invariant quantities in equations 
(I26I31|) and H36I38|I . but constructed as a metric and energy-momentum perturbation of the pulsating star. Then, we 
can expand in A to get 



(83) 



It is important to remark here that the (1,1) superscript refers not only to quantities constructed from the g^^'^) 
perturbations, but in general to any perturbative quantity of order Ae. In general Q^^'^'^ and Q^^'^'^ can be expressed 
as: 



^(0,1) ^ -^(0,1) ^ 

^(1.1) ^ 7^(1,1) +y 1(1,0) 



Jr 



(0,1) 



(84) 
(85) 



where the objects H^^'^'' and are linear in the (0, 1) perturbations, while ji^'^^ and H^^'^"^ are respectively linear 

in the (1,0) and the (1,1) variables. It is clear that 7Y(0'1) is nothing but any of the gauge-invariant quantities in 
(|26I31|) and (|36I38I) for the special case of a static background, as it must be for the (0, 1) perturbations. The quantity 
7-^(14) jg constructed from (1, 1) quantities in the same way but, as we are going to see, it is not a gauge-invariant 
object at the order (1, 1): we have to add some extra terms of the form given in H85|) and which come from the ansatz 
(|83|l . In what follows we show that this gives the desired gauge-invariant (1, 1) quantities. 

Gauge transformations and gauge invariance in 2-parameter perturbation theory have been studied in |l7l IT^ . The 
gauge transformation of a first order perturbation of order (0, 1) of a generic tensorial quantity T is given, as in H25|l. 
by 



£c r("^°) 

S(o,i) 



^(0,1) ^ ^(0,1 

whereas a second order perturbation of T, in particular of order (1,1), transforms according to 

t(M) ^ + £^^^^^^ T^m + T(^A) + (^^^^^^^ + , ^^^^^ J) r(o.o) 



(86) 



where { , } stands for the anti-commutator {a, b} = ab + ba. 

In the present case, we have chosen to fix the gauge for the radial perturbations g'^'^-*, 
simplifies the previous transformation rule to 

^ ~ ^ ^ 4(0,1) ^ ^ 
We have assumed that C/("^^Hs a gauge-invariant quantity at order (0, 1), therefore 

^(0,1) ^(0,1) ^ -^(0,1) _ -^(04) ^ . 

From (I85|l and the fact that we have fixed the gauge for radial perturbations, we can write 



(87) 

see equation (|Ffj) . This 



a 

Furthermore, we note that every H(i'i) and can be expressed as follows: 



-J. 



(0,1) 



(1,1) 



(0,1) 



.(0,1) 



(89) 



(90) 



(91) 



where A and are linear operators involving differentiation with respect to the coordinates of Af ^ and integration 
on 5^. These operators act on spacetime objects and return objects with indices on M^. From now on, for the 
sake of simplicity, we only consider metric perturbations. The corresponding procedure for energy-momentum tensor 
perturbations follows the same lines and is given in Appendix [CI Using the gauge transformations 1)86(1 and H88|l . the 
transformation rules for 7i(^'^) and are given by 



7^(1,1) 

J(0,1) 



7^(1,1) 

J(0,1) 



A 
B, 



£ 



r(l,0) 



«(0. 

4(0, i)& 



(92) 
(93) 
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Moreover, we know pit Is^ that the quantities in (|26I31|) are gauge-invariant as first order perturbations, hence 
A [£^g^^'^^] must vanish for any vector field f . Therefore, (|92|l becomes 



°^?(0,l)O 



7:^(14) ^^(1,1)^^ 

and the gauge transformation reduces to 

G^'-'' - +^ r^^^^^^gd^o)! +^1(1,0)^^ kc.g^"''^^ 



(94) 



(95) 



Using this expression we can show that the variables C/*-^'^-' are gauge-invariant. For the sake of brevity we only 



give here the proof for the metric perturbation fc^^ H2t)|) , the analysis for the other metric and fluid perturbations 
m mm and mm is given in Appendix El To proceed with the proof we expand the generator of the gauge 
transformations associated with the non-radial perturbations in tensor harmonics. 



(96) 



Since the metric perturbations g^^'") do not depend on the coordinates of S^, and taking into account the gauge choice 
(|57|l . we have 



■^5(0,1) &Aa — "'AC ? ^ ■'^ 



£ 



5(0,1) 6ah 



0. 



(97) 



where £^ is the Lie derivative acting on M^. Now, we apply the ansatz described in to Hab and we get the 
following expressions 



^(1,1) _ . (1,1) 

''■AB — "-AB 



(1,1) , (1,1) 

Pa\b 



^(1,0)C _ „p(l,0)C 
■^AB — '^'^ AB ' 



7(0,1) _ (0,1) 



(98) 



where pA is defined in 1)321) and F^^'*'^ are the radial perturbations of the ChristofFcl symbols. From (|94() and the 



analysis of |21|, |33 , we have 



-£/(l,l) _ ay(l,l) , /. 
''■AB — ''■AB ^-^^ 



1,0) _ / (1,0) (1,0) 
Ms \Pa\B ^ l^B\A 



where the explicit expression of the Lie derivatives is 



£ih 



(1,0) 



^"'AB 



c,,(i,o) 

AB\C 



;,(i,o);c 

"CB ^\A 



,(1,0) £C 
^CA 



C|B I ^^Pa ^ "'AC ^ 



,(l,0)eC 



(99) 



(100) 



Finally, introducing all the expressions (|98I100|I into the gauge transformation law 195|l . we get the gauge invariance 
of kAB at (1, 1) order: 



r(i,i) _ , (1,1) 

'^AB — '^AB 



(101) 



C. Equations for the Ae coupling perturbations 

The explicit form of the equations that govern the behaviour of the coupling terms is obtained by introducing in 
equations (|42I50() the following expressions: i) for the background quantities we will use the expressions of the GSGM 
quantities describing the radially oscillating spacetime (the static background plus radial perturbations), given by 
equations H75I82|I : ii) for the perturbative quantities we use the corrections to the radially oscillating star, that is, the 

quantities that come from perturbative terms like g^^^^ = gl^"^^^ -I- -^g^^^^ ■ Once we have introduced all these quantities, 
expanded the equations and extracted the Ae part, we get a set of equations that can be expressed as a linear non- 
radial operator Z/nr acting on the (1, 1) variables, and a source term <S built from the (1,0), (0, 1) quantities, see 
Eq. ®. 

As we explained in Section m this particular structure of the (1,1) equations is quite convenient in order to build 
an initial-boundary value problem and solve it numerically by using time-domain methods. The basic idea is that 
given a numerical algorithm capable of evolving linear non-radial perturbations, we can build an algorithm for our 
(1, 1) perturbations by just adding source terms to the original algorithm. 
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The time evolution of non-radial perturbations of a static star has been successfully analyzed by numerically 
integrating different systems of perturbation equations [s^ |3 . Taking into account the main features of our 
formulation, the scheme introduced in ^ seems to be adequate for the purpose of implementing a numerical code 
to solve our perturbation equations. One of the main points in the scheme introduced in |39l l4lj is the fact that the 
Hamiltonian constraint is not just an error estimator for the evolution equations, as it is usually done in many free 
evolution schemes. In the scheme of |39t l4lj |. the Hamiltonian constraint is part of the system of equations and it is 
solved at every time step for the perturbative quantity k, Eq. H27|) . This provides some control of the errors induced 
by constraint violation. As a consequence, the resulting numerical code 39] is able to evolve non-radial perturbations 
for long times and is capable to estimate the damping time and mode frequencies with an accuracy comparable to 
frequency domain calculations. 

The main idea of our present ongoing work 10] on the numerical solution of our perturbative equations is to follow 
the scheme of [s^ ^] . Taking into account our discussion in Section ^1 and above about the general structure of 
the perturbative equations, in particular the same differential structure of the perturbative equations at the orders 
(0, 1) and (1, 1), it is clear that this scheme is easily portable to our problem. To that end, it is very important the 
fact that the Hamiltonian constraint is solved for one of the perturbative quantities since at every time step we need 
to evolve the equations for the (0, 1) and (1, 1) perturbations. This means that if we do not solve the Hamiltonian 
constraint, the errors accumulated from constraint violation would be double than in a standard computation of 
non-radial perturbations. Therefore, the use of the scheme of j3^ is a key ingredient in trying to obtain accurate 
long term evolutions. We expect that the resulting numerical code would allow us to investigate the non- linear effects 
of coupling. In particular, we are interested in looking for non-linear harmonics, possible resonances, parameter 
amplification, and/or changes in the damping time of non-radial perturbations. 

In the stellar interior, we evolve the (hyperbolic) equations for the metric perturbation -x^^^'^") and for the fluid 
perturbation H^-^'^\ which in some particular gauges coincides with the enthalpy perturbation. The Hamiltonian 
constraint provides us the metric perturbation k^^'^\ Subsequently, all the other metric ^^■^'^"> and fluid (7'"'^'^^ a^^^^-*) 
perturbations can be obtained from the perturbative equations H44|l . H49|l and (|45|l . 

The wave equation for x'^^''^'' and the Hamiltonian constraint are given by (|42l) and USUI) respectively, while the 
sound wave equation for H^^'^'> has to be determined. We define the fiuid perturbation ff^^^ (sec Appendix [21 for a 
proof of the gauge invariant character of this quantity) as 



2(0) (0) 

s P (1) 



(102) 



p(") + pi") 

where the superscripts (0) and (1) have the meaning already explained in Section^ The sound speed in the radially 
pulsating spacetime can be split as follows 



dp 



(103) 



In particular gauges, the Regge- Wheeler ^3] one for instance, the gauge-invariant quantity uj^^^ coincides with the 
gauge dependent perturbation uji^) [see Eq. H38ll ]. and Hi^) describes the enthalpy perturbation. 



^(1) 



p(") + 



(104) 



where 5pi^) is defined by H35|l . The wave equation for H^-^'^) is obtained as a linear combination of the time frame 
derivative of equation H47(l and the spatial frame derivative of (|48|l . After having introduced the equations l|44l 1461 
1481 1431 150|l to reduce the number of perturbative unknowns and the transformation H102|l , we have the following wave 
equation (written in the GSGM formalism): 



- H 



clH" 



^ff = , 



(105) 



where Th contains all the remaining terms (with derivatives of lower order) . The complete equation has been written 
in Appendix^ The wave equation H105|l is valid in the GSGM framework for barotropic non-radial perturbations on 
a time dependent background. In case of a static background, provided the introduction of the background quantities 
(|55l56f) . it reduces to an equation well known in the literature (see i.e. j3^, (3^, [39|). 
We can now write the perturbative equations for the stellar interior. We consider instead of the perturbative 



quantity which diverges like r as we approach spatial infinity, the perturbation variable S*'^'^) 

of course is well behaved at infinity. This quantity satisfies the following gravitational wave equation: 



,(1,1) 



/r which 



0(1,1) 2(<I.-A) 0(1,1) 

'~',tt I ^ '-',rr 



.2(*-A) 



4 /I 



^2A 



14 



- (5 + 4$,rr) + 



(/(/ + l) + 2)e2A 



14) 



2* c 

e 6s . 



(106) 



where Ss denotes the source term for this wave equation. The source terms in our (1,1) perturbative equations have 
the following pattern 



(107) 



which shows how the source terms introduce the coupling between radial and non-radial perturbations in the (1,1) 
equations. In particular, the source term in the gravitational-wave equation, Ss has the following form 



0-5 



2e 



where the coefficients at are just linear combinations of radial perturbations with coefficients constructed from back- 
ground quantities. Their explicit form is given in Appendix IXI 

The perturbative fluid variable H^^'^^ also satisfies a wave equation, but with a different propagation speed. We 
call this equation the sound wave equation. It has the following form: 



(l + 3cl) {A^r + 



2 

- + 2$, 
r 



c, e 



-A,, 
^(1.1) 



1.(1,1) 



-2$2^ + [(3$,, + A,,) r + 1 - e^A] ^ 



H 



(109) 



and the source term can written as 

r(0,l) , I, u(0,l 



,(0,1) 



+ b, (r5(°-i) + fc("'i)) + 5io7.l"''^ + &ii7^°''^ + b,2^l^f^^ + 613^^°^'^ + Wac 



(^5(0,1; 



,(0,1) 



(110) 



where the coefficients bi have the same structure are the coefficients in (|108|) . Their explicit expressions can be 
found in the Appendix 1X1 

For the last perturbative variable, the metric perturbation k^^'^'> , we will use the Hamiltonian constraint instead 
of an evolution equation. After some calculations we get: 



M) _ ^(;,i) + ("2 _ ^ \ fc(i,i) _,. 2 ^^^^ ^ ^ ^(1,1) _^ 1 _ ^ g2A ^ 2A ,r - 1] fc(i'i) 



1 [1(^1 + i)e2A + 4 - 4A,,r] ^^-D = Sna. 



(Ill) 



where Snamii is the source term for the Hamiltonian constraint. As in the previous equations, it follows the pat- 
tern p07|) . The precise form of SHamii is: 



(112) 



The coefficients Ci, in the same way as the coefficients and bi only contain radial perturbations 5(1^'^) and quantities 
associated with the static background. They are also given in Appendix ^ It is worth to remark that the polar 
non-radial perturbation equations on a static background are obtained from equations (|106I109I111|) by discarding 
the source terms and replacing all the (1, 1) perturbations with the corresponding non-radial (0, 1). The sources are 
determined from first order perturbations. The radial perturbations from the equations H64I68|I . and the non-radial 
perturbations (described by the quantities S"*-"'^-* , k^'^'^^ , and if^°'"'^') from the first order analogous of the above system 
(see ll^), and the equations if^B^HS)) adapted to a static background to get the ^/;(°'^\ j(o,i} g^^^ a^'^'^K 

The stellar exterior is described by a Schwarzschild spacetime on which gravitational waves carry away some energy 
of the stellar oscillations. All fluid perturbations are not defined outside the star and the radial perturbations vanish 
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because of Birkhoff 's theorem. Therefore, the source terms in our perturbation equations vanish. Only the metric 
perturbations survive, and they satisfy the gravitational wave equation (|106|l and the Hamiltonian constraint (|lll|l . 
which take the following form: 



o2(<E'- 



6M 



1,1) 



2M 



2M 



,2A 



l{l + l) 



^(14) 



4M 



M 



,2A 



0, 

(113) 



-2A 



114) 



It is worth to mention that the above equations coincide with the equations for non-radial perturbations of a static 
stellar background outside the star, as expected. 

On the other hand, Zerilli showed that the even-parity perturbations of a Schwarzschild background have just 
one degree of freedom, and therefore can be described by just one variable, the Zerilli function, satisfying a wave 
equation. At order (1, 1) the Zerilli function can be built from the two metric perturbations S"'^'^^ and fc^^'""^' and their 
derivatives, as at first order ^^nd is given by 



^(1,1) 



2r^e 



{l + 2){l-l)r + 6M 
It satisfies the Zerilli equation 0, 

- + e2(*-^)ziV) -I- ^e^^Zt^'i) _ v-(^)z(i,i) = q 

where V{r) is the Zerilli potential [45l |: 

V{r) = 



2M\ ni{ni ~ 2fr^ + 6(n, - 2fMr'^ + 36(n, - 2)M^r + 72M^ 
r^[{ni - 2)r + 6M]^ 



r / 



(115) 



(116) 



(117) 



where the quantity ni ~ 1(1 + I) has been introduced to simplify the expression. 

Finally, we can determine the power of the gravitational radiation emission at infinity by using the following 
expression 



dE 

H 64tt 



1 ^ [1 + 2)1 ^ ,2 



I. , 



il^2)\ 



(118) 



VI. BOUNDARY CONDITIONS 



In this Section we discuss the boundary conditions at the origin and at the stellar surface for the (1, 1) perturbations 
describing the coupling of radial and non-radial modes. With regard to the outer boundary, by locating it far enough 
we can use the well-known Sommerfeld outgoing boundary conditions on our fields. 

At the origin, the boundary conditions are just regularity conditions on the perturbative fields, which can be 
obtained by a careful analysis of the equations that they satisfy. The analysis of Taylor expansions of the differential 
operators that a ppe ar in our equations near the origin leads to the following behaviour for the non-radial perturbations 
g(o,i) and g(i'i) M- 

I > 0, 5(i)^r'+i V^i)~r'+i, (119) 

I > 1, 7(1) ~ r'-i i/^i) - r' a-r'. (120) 

By using the behaviour given by these expressions and also the behaviour given in expressions (|72l) for the radial 
perturbations it can be proved the regularity of the source terms at the origin. 

At the surface of the star, we have to consider a boundary condition for the matter variable H, which vanishes in 
the spacetime region outside the star. For the metric perturbation variables S and A:, since they must be continuous 
through the stellar surface, we can just use the junction conditions at the surface to determine them. 
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To that end, let E be the surface of the static unperturbed star (i.e. r = Rs). The surface of the perturbed star 
can then be described in the foUowing way 

E = {.T + A^(i^o) + e^("^i) + AeC^i'i) : a; e e} , (121) 

where S}-^ '■'^ is a vector field that denotes the Lagrangian displacement of a fluid element due to the action of per- 
turbations of the order {i , j). A physical requirement that follows from junction conditions is the vanishing of the 
unperturbed pressure p at the unperturbed surface E. In the same way, the corresponding boundary condition for the 
perturbed spacetime is the vanishing of the total pressure p + AJp*-^'"' + e<5p'-°'^' + Ae^p'"^'^^ at the perturbed surface 
E. This condition turns out to be equivalent to the vanishing of the Lagrangian pressure perturbations on E, the 
unperturbed surface, at every order. The Lagrangian pressure perturbations are given by: 

Ap(i'O) = <5p(i'0) + £^,,„,p, (122) 
Ap(M) ^ ^p(o.i) + £^^^^^p, (123) 

= + ~ i , }^ p, (124) 

where 5 and A denote the Eulerian and Lagrangian perturbations respectively, and we have used the lower order 
boundary conditions Ap(^'°^ = Ap(°'^^ = in order to simplify the condition p24|l . 

From this analysis we can conclude that the boundary conditions for the fluid perturbations are described by the 
set of expressions given in p22ll24|) . However, in practice, in many applications of first order perturbation theory, 
dynamical boundary conditions cither for density or enthalpy perturbations have been considered. This alternative 
boundary conditions follow from the analysis of the time derivative of the condition p23|l (see [s^l for more details). 
In our current development of the numerical implementation of the perturbative equations we are considering both 
types of boundary conditions with the perspective of analyzing which type works best for our formulation. 

Finally, the junction conditions for the metric perturbations can be determined by i mpo sing continuity of first 
and second fundamental differential forms and their perturbations at the surface |47l liM liil l5(l | . The explicit form 
of these conditions has been presented in |23| for first order perturbations of a time-dependent stellar background. 
Alternatively, one may use the "extraction formulas" pol | that relate the Zerilli function with metric perturbations at 
the stellar boundary. 



VII. REMARKS AND CONCLUSIONS 



Non linearity is the rule rather than the exception in dynamical phenomena in all branches of physics. The 
modeling of compact objects such as neutron stars and supernova; core must ultimately be rooted in general relativity 
(or some of its generalisations), where non-linearity represents a fundamental physical character of the theory, i.e. the 
self-interaction of the gravitational field, and not just corrections to an underlying linear modeling of gravitational 
phenomena. In relativistic theories of gravity, gravitational radiation is the typical outcome of dynamical phases in 
the life of sources such as binary systems and supernovse, and major experimental efforts are currently under way to 
detect this most elusive prediction of Einstein gravity for the first time. This will eventually lead to the development 
of a whole new branch of astronomy, based on observing gravitational radiation, much in the same way it has been 
in the past for x-rays and other parts of the electromagnetic spectrum outside the visible band. In this context the 
accurate theoretical modeling of sources is crucial to the final end of providing templates in this game of looking for a 
needle - the signal - in the haystack, the noise. While ultimately a full numerical relativity description of gravitational 
wave sources is needed to model the most non-linear dynamical phases, much interesting physics can be understood 
by using approximate methods. Furthermore, a semi-analytical approach typically helps to shed light on the physical 
processes, thus complementing the numerical work. 

Relativistic perturbation theory is ideal for those cases where a known solution of the field equations is explicitly 
known, as for black holes, or can easily be obtained, as is the case for compact stars. An advantage of the relativistic 
perturbative approach is that it directly incorporates gravitational waves. For smaller perturbations, linear theory 
suffices. If one wants to consider mildly non-linear oscillations of a compact object, second order effects and mode 
coupling have to be taken into account. For black holes, many studies already exists in this direction (see e.g. 

llSJi^). In the case of neutron stars, while linear perturbations and instabilities have been studied for long 
time m, relatively little is known of non-linear dynamical effects, mostly through numerical studies (see e.g. 
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13. Isll^. ItI Is^ ) . A second order perturbative approach is therefore timely and may help to understand known 
problems and even reveal a new phenomenology. 

In this paper we have developed the relativistic formalism to study a particular second order effect, the coupling 
of radial and non-radial first order perturbations of a compact spherical star. From a mathematical point of view it 
is very convenient to treat the two sets of perturbations, radial and non-radial, as separately parametrized, using the 
multi-parameter perturbative formalism developed in |l7L Il8l|. Then we have considered the expansion of the metric, 
the energy-momentum tensor and Einstein equations in terms of two parameters A and e, where A parametrizes the 
radial modes, e the non-radial perturbations, and the Ae terms describe the coupling. This approach provides a 
well-defined framework to consider the gauge dependence of perturbations. In this mathematical context we have 
imported the formalism of Gundlach and Martin Garcia 0, and Gerlach and Sengupta |2ll , describing gauge- 
invariant perturbations of a general time-dependent spherical spacctime, expanding the latter in a static background 
and a radial perturbation. Fixing the gauge for radial perturbations allows us to: i) use the GSGM gauge- invariant 
non-radial e variables on the static background; ii) define new second order Ae variables, describing the non-linear 
coupling of the the radial and non-radial linear perturbations, that are also gauge-invariant at the Ae second order. 
This higher order gauge invariance, attained by partially fixing the gauge at first order, is similar to that considered for 
example in (22] and ^3]. In our case however we use a 2-parameter A — e expansion |l7lll8l|. so that we only need to fix 
the gauge for radial perturbations. Assuming a barotropic perfect fluid, we have derived the evolution and constraint 
equations for our variables, in particular those for the coupling terms of order Ae, focusing on polar perturbations. We 
leave for future studies the implementation of more realistic equations of state, such as the non-isentropic one used in 
[5^ Is^ l . As expected, in the interior the Ae variables satisfy inhomogeneous linear equations where the homogeneous 
part is governed by the same linear operator acting on the first order e non-radial perturbations, while the source 
terms are quadratic and made of products of A and e terms. In the exterior there is no direct coupling, and the whole 
dynamics is embodied in the Ae order Zerilli function. Thus the effect of the coupling is transmitted from the interior 
to the exterior through the junction conditions at the surface of the star. Finally, we have given a brief discussion of 
the boundary conditions, focusing on those on the surface. These are typically expressed in terms of the Lagrangian 
pressure perturbation, therefore we have defined a Ae second order Lagrangian displacement and a corresponding Ae 
Lagrangian pressure perturbation, appropriately related to the Eulerian perturbation. Thanks to the vanishing of the 
first order A and e Lagrangian pressure perturbations on the surface, this relation turns out to be linear. 

Work is currently under way, numerically implementing the formalism presented here, in order to provide a first 
analysis of the possible effects of the coupling between radial and polar non-radial perturbations ^10] . Some of these 
effects are easily anticipated for the case of axial oscillations. These are decoupled from fluid perturbations at first 
order, but are driven by the radial pulsations at the Ae order |^. Eventually these studies may possibly even lead 
to discover new unexpected effects of mode coupling. Surely we expect to find non-linear harmonics arising from the 
radial non-radial coupling, similar to those between various radial modes found in [s^ for tori around black holes, 
a prediction that appears to be confirmed by a numerical relativity study of neutron stars in the conformally-flat 
spacetime approximation [ssf . 
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APPENDIX A: SOURCE TERMS FOR THE (1, 1) PERTURBATION EQUATIONS 

In this Section we give the expressions of the coefficients appearing in the source terms of the equations (|106I112|I for 
the (1,1) perturbations. 

In the case of the gravitational wave equation (|108|l the source term has the form: 

Ss - ai5,i"'i) + a^S^^^'^ + aa^,^'^ + a,S^°''^ + a, ~ 2e^-*fc(°'^)) + a,k^"^'^ + ar^^^'^^ , (Al) 

where the coefficients are given by 

ai=2(r5(i^")-77(i^"))e-2'^, (A2) 
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03 
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For the sound wave equation the source term is given by Hll()(l : 



/,(0,1) 



and the expression of the coefficients hi is the fohowing 
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5i4 = -- 



p r 
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2/(; + l)7(i'")e*}, 
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Finally, in the case of the Hamiltonian constraint, the equation we have considered for fc^^'^\ the source term is 
given by (|112|l : 



where the coefficients Ci are given by 
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APPENDIX B: SOUND WAVE EQUATION 



We write here the complete sound wave equation p05|) of a generic barotropic fluid perturbation on a time-dependent 
background in terms of the GSGM quantities, 
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APPENDIX C: GAUGE IN VARIANCE 



In this Section we show the gauge-invariant character of the perturbative quantities (|27I31() and H36I38|I at the 
perturbative order (1, 1). To that end we foUow the procedure described in Section I V ij| Namely, we determine, for 

each quantity, the corresponding term 'Y^^x'^''^^ j'^''^^ in the expansion (|85|) and the gauge transformation of liS^'^^ 
[see Eq. (|94|) ]. Then, considering the gauge transformations for the non-radial perturbations J^^"'^-' and introducing 
all the terms in (|95|l . we will prove the gauge invariance of our perturbations. 
The non-radial gauge transformations (which we will use later) are 

f/^ = vT^ + Ia , G("'i) = G("'i) + 2e. (CI) 

Metric perturbations. In Section ^ we have shown the gauge invariance of the metric perturbation k^^g"^ ■ Here we 
prove that of the metric perturbation k^^'^\ From H97I) we find that, 

^^(1,0)^(0,1) ^2g(i,o)AB^^p(o.i)^ (C2) 
^(1,1) ^ ^(1,1) „ 2v^g'i;^k^' , (C3) 
where va = f\A/''' ■ Therefore using (inil) we find from ^ that fc(i^i) = fc^^'i) . 

Stress-energy tensor perturbations. They are described by the seven quantities (|28I31() . The corresponding terms 
E.I^'^'^J^"^'^ are: 



: + 2Qva) g^'^^^^^pf'^ - (g[/) + 2Q(i^°).^) g^^pf'^ + fcll q(i.o)g(04) , (C5) 

Ta : (4''"^'' - tABg^'''^^^) p'^''' - yQ(^'°)G[°'^) , (C6) 
: r2Q(i^o)G(°'i) . (C7) 

where Q is the pressure. The gauge transformation of the W^^'^^ part of the quantities under discussion is given by 
^i(o,i)^^af^ ^ where the energy-momentum tensor for radial perturbations has a block diagonal form 



Then, we find 



&°)=diag(4^f ;r2Q(i^o)7,,) . (C8) 



£i,o.jAf = (i^^c°^l^ + r^Q(^-°^e|A)Fa, (CIO) 
£i,oJaf^ = iQ^^'^i'' " Kl + + 2vcQ^''"^i^) r7a. + {2r'Q^'''k) Z., . (Cll) 
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Therefore the gauge transformations of the Ti.^^'^^ terms for the above quantities are: 

(C12) 

QvA)g%^'^^i\ (C13) 

(CM) 
(C15) 

Combining all this terms into the relation (|95|l . we finish the proof of the gauge invariance of (|28I31I) at the 
perturbative order (1,1). 

The fluid perturbations for a barotropic fluid that have to be considered are the two components of the velocity 
(|36l37f) and the energy density (|38|) . Applying the same procedure we find the corresponding J^}''^'' terms 
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The gauge transformations (|94|l for Ti'^'^^ are now determined by the Lie derivative of the radial perturbations of the 
fluid velocity and energy density 



where 57 = In p. We obtain 
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The corresponding liP-'^^ transformations are 
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It is again easy to verify the gauge invariance of these perturbations by bringing all terms into the relation (|95|l . The 
fluid perturbation H^^\ defined in Eq. (|102|) . can be expanded like the previous quantities, 
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Therefore, the gauge-invariant character of and H^^'^\ having fixed the gauge for radial perturbations, follows 

from the gauge invariance of w^'^'^^ and which has already been proved previously. 

APPENDIX D: CONNECTION TO THE REGGE- WHEELER METRIC VARIABLES 

The perturbations of spherical stars and Black Holes are commonly studied in the Regge- Wheeler (RW) gauge 01 ■ 
Therefore, in this Appendix we provide the relations between the perturbative variables used in this work and those 
of RW. In the RW gauge the linear polar non-radial perturbations assume the following expansion in spherical tensor 
harmonics, 
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where i?o°'^, Hf'^, ^2*'im' ^im^^ ^'^^ functions of {t,r). The Einstein equations for a spherical star imply that 
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im ~ ^-"2 im- choose the RW gauge for the perturbative variables describing the coupling by enforcing 

the RW form of the polar metric perturbations at the (1, 1) perturbative order. Then, we impose the following form 
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where also i?!^;!^ i?/^'^^ 



o,im' --i,im' ^^2.im^ ^hn^' are fuuctious of {t,r). 



The expansion of the linear non-radial and the coupling metric perturbations in terms of the GSGM variables can 
be derived from (|23|) by applying the 2-parameter expansion of GSGM formalism described in section^ First of all, 
we must first apply the RW gauge at first and second perturbation order, i.e. 
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Then, we take the perturbative expansion of the gauge- invariant tensor Uab I|39(I , written in the basis of spanned 
by the vectors H75l7t)|) . and use the definition 14U|I for the x perturbation. Finally, taking into account Einstein's 
equation 1)41(1 . we find the following relations: (i) For the (0, 1) linear perturbations. 
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(ii) For the (1,1) coupling perturbations, 
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Before concluding this section we also give here the relations even between the quantities used in this paper and those 
considered by Allen et al. and Ruoff [3| at linear order: 

(0,1) _ / S' in [33] , (o,i) _ ( F/r in \^ , , . 

^ ~\rS in [3^, " -\r/r in[34]. ^ ' 

For more details about the connection between the GSGM and other formalisms see the appendix in [l^ . 
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